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Abstract 

We study a two-loop induced radiative neutrino model at TeV scale with global U{1) symmetry, in which 
we analyze dark matter and resonant leptogenesis. The model includes two kinds of dark matter candidates. 
We discuss what kind of dark matter can satisfy the observed relic density as well as the current direct 
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resonant leptogenesis can be differentiated from the other scenarios at TeV scale or not. 
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I. INTRODUCTION 


After the discovery of the Higgs boson at the LHC, the Standard Model (SM) has been estab¬ 
lished well. However the SM still has to be extended in order to explain the existence of dark 
matter (DM), the small neutrino masses, the baryon asymmetry in the universe and so on. Radia¬ 
tive seesaw scenarios are renowned as one of the economical models which simultaneously explain 
the existence of DM and the neutrino masses. Since the DM candidate is necessary to generate the 
neutrino masses in this kind of models, physics between DM and neutrinos is strongly correlated in 
this simple framework. For example, couplings and mass scale of DM are related with the neutrino 
mass scale. In addition, since this kind of models can naturally include a new particle with TeV 
scale mass, radiative seesaw scenarios have good testability in near future experiments. Along this 
line of idea, a vast literature has recently arisen in Ref. [1-72]. 

On the other hand, explaining the observed baryon asymmetry in the universe via leptogenesis 
is one of the challenging issues in the framework of radiative seesaw models with right-handed 
neutrinos, since couplings related to the source of the leptogenesis are expected to be 0(1) due 
to the requirement of the neutrino masses. It causes the strong washout of the generated baryon 
asymmetry. In order to avoid this problem, we have to take hierarchical couplings with highly 
degenerated masses between the source and the mediated fields.^ For example, generating the 
baryon asymmetry via resonant leptogenesis has been discussed based on the Ma model [73]. 

In this paper, we study a two-loop induced radiative neutrino model at TeV scale with a global 
U{1) symmetry, in which we analyze DM and resonant leptogenesis simultaneously. In this model, 
we have a scalar or a fermion DM candidate. We discuss which kind of DM candidate can satisfy 
the observed relic density as well as the current direct detection bound, and can also be compatible 
with leptogenesis. Since our model has two sources of leptogenesis, we also show different points 
of our resonant leptogenesis from the other scenarios such as the Ma model at TeV scale [73]. 

This paper is organized as follows. In Sec. H, we show our model including field contents 
and their global U(l) charges, Higgs potential, and neutrino masses. In Sec. HI, DM properties 
including relic density and current limit by direct detection experiments are discussed. In Sec. IV, 
we analyze resonant leptogenesis. Summary and conclusions are given in Sec. V. 

^ Such couplings can be achieved by making use of the experimental fact that one of three neutrino masses can be 
negligible. To get the sufficient baryon asymmetry via thermal leptogenesis in the radiative seesaw framework, 
resonant leptogenesis would be only the solution that requires the mass degeneracy between the source and the 
mediated fields. 
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II. THE MODEL 


A. Model setup 
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TABLE I: Field contents and charge assignments of SU{2)l x U{1)y x C/(1), where indices i = 1 — 3 and 
j = 1,2, (3) represent the generation. 

As shown in Tab. I, we introduce two (or three) gauge singlet vector-like fermions and 

gauge singlet Majorana fermions Nji, and Xr as new fermions. The number of these particles 
should be more than two in order to obtain at least two non-zero neutrino mass eigenvalues. We 
also introduce an inert SU{2)l doublet scalar rj, two neutral inert singlet scalars {x^^x^')i a 
neutral singlet scalar E as new scalars. We assume that only the Higgs doublet field in the SM $ 
and the new singlet scalar E have vacuum expectation values (VEVs), which are symbolized by 
($) = u/\/2 and (E) = v' jy/2 respectively.^ We impose a global U{1) symmetry, under which 
does not have a charge in order not to couple to the Goldstone boson (GB) [37]. The global ?7(1) 
charge x 7 ^ 0 is in principle an arbitrary, and the field assignments play a crucial role in realizing 
our neutrino masses at two-loop level. If the 17(1) charge x is fixed to be x = 2, we can identify 
this 17(1) symmetry as a global B — L symmetry. Hereafter we assume this global 17(1) symmetry 
to be a kind of U{\)b-l symmetry. Note that while the new fermions are added as vector-like 
and do not contribute to anomalies, this model is anomalous since the three chiral fermions with 
B — L = —1 corresponding to the right-handed neutrinos are not introduced. If one would like 
to have an anomaly free model, the anomalies can be cancelled by introducing some pairs of new 
heavy vector-like fermions [74, 75]. However this is beyond the scope of this paper. This model has 
an accidental Z 2 symmetry which can assure the DM stability, and the accidental Z 2 assignments 
are shown in Tab. I. 

The renormalizable Lagrangian for Yukawa sector, mass term, and scalar potential under the 

^ The scale of v' should be larger than v' ~ 10^ GeV for successful leptogenesis as we will see later. Otherwise 
the annihilation channel NiNi, XiXi —> GG whose reaction rate is determined by v' does not satisfy the out-of¬ 
equilibrium condition at T ~ 0(10) TeV where T is the temperature of the universe. Thus the baryon asymmetry 
would be washed out due to this process. 
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where the first term in Cy generates the SM charged lepton masses, and we assume the couplings A, 
A', \" and in the scalar potential to be real for simplicity. As we will see below, these couplings 
become important for neutrino mass generation. 

After the symmetry breaking, the scalar fields can be parametrized by 


^=( V^=f V^ = (II.3) 

\j^iv + <P + iz)J \^iVR + iVi)J v2 

where v ~ 246 GeV is the VEV of the SM Higgs doublet, and and z are respectively the GBs 
which are absorbed by the longitudinal components of the W and Z bosons. Inserting the tadpole 
conditions, the resulting mass matrix of the GP even scalar (cj), a) is given by 


rr?{(t),a) = 


A$u^ 2\^y,vv' 
2\^yvv' Xyv''^ 


cos a sm a 
— sin a cos a 


mi 0 

0 m'ij 
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(11.4) 


where h is the SM-like Higgs boson and H is an additional CP-even Higgs mass eigenstate. The 
gauge eigenstates (j) and a are rewritten in terms of the mass eigenstates h and H as 



with sin 2a = 


4A<j)SUti^ 
m]j-ml 


(H.5) 


The GB G in Eq. (H.3) appears due to the spontaneous symmetry breaking of the global 17(1) 
symmetry. The couplings between the GB and the particles with non-trivial global 17(1) charges are 
given by J^d^G/v' through the global 17(1) current J^. As one can see, the coupling is suppressed 
by the VEV v'. 

The mass matrices of the CP even and CP odd states of the inert scalar bosons {y,X^■:X^') r/i 
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are respectively given by 
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where we define diagonal mass matrices and their mixing 

matrices so that they satisfy . Depending on the couplings, 

the lightest CP even or CP odd mass eigenstate with the mass mn^ or mi^ can be a DM candidate.^ 
The non-standard couplings between DM and the GB induced by the non-self-conjugate couplings 
A, A', X" and may be relevant to compute the DM relic density. This coupling can be written 
down as 


I ^''"r^Rr^R , r^R/^R 


VdM-DM-G-G — — ( -^—,Oii 02 i + — —O11O31 -I- —O2W31 + 
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(II. 8 ) 


' 4\/2u 

with the mixing matrix O^. In addition, the couplings between the CP even Higgs bosons and 
the GB are also relevant to compute the DM relic density. These couplings come from the kinetic 
term of S and can be written as 


C D 


— sin ah + cos aH {—sin ah + cos aH) 
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2 u '2 
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(II.9) 


Finally the mass eigenvalue of the charged inert scalar 77 + is given by 

A$^u^ -I- Xr^Y.v''^ 


2 2 , 

m^± = m^ + 


(II.IO) 


In this model, the typical mass scale of these new exotic particles is assumed to be TeV scale. On 
the other hand, we should take v' ~ 10^ GeV for successful leptogenesis. Therefore a certain degree 
of tuning among the relevant couplings cannot be avoided. More specifically, demanding that the 
diagonal elements of the mass matrix Eq. (II.4), (II. 6 ) and (II.7) are TeV scale and off-diagonal 
elements are 10 GeV scale to obtain small mixings of the order of ~ 10“^ {i ^ j), the order 
of magnitude of the couplings should roughly be X" ~ 10“^^ A, A', A$s ~ 10“’^ and A^s, A^s, 
~ 10“®. Although this point may be a disadvantage of this model, it would be worth 


^ As we will discuss later, the CP even state is identified as DM. 
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discussing such a new concrete model with a global U{1) symmetry as an example model since all 
the phenomenology of the neutrino masses, the existence of DM and the baryon asymmetry of the 
universe are closely correlated. 


B. Neutrino mass matrix 


Due to renormalizability and the strong restriction of interactions via the global U{\) symmetry 
in this model, neutrino masses are not generated neither tree level nor one-loop level. If the vector 
like fermion F has a Majorana mass term, neutrino masses would be generated at one-loop level 
(for example see Ref. [1]), however this is not our case. As a result, neutrino masses are induced at 
two-loop level, and we have three types of diagrams as shown in Fig. 1. The formula of the total 
neutrino mass matrix can be given by 
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where a pair of and a pair of ^ and are the dimensionless 
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loop functions which come from the left, center and right diagrams in Fig. 1 respectively. These 
loop functions are defined by 
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FIG. 1: Radiative generation of neutrino masses. 


where 


nR nR (1^ 

am'^ bm cn'^ dn 


and 


jiijk) 

l(Rmtn) 


jiijk) 

2R{RmIu) 


riijk) 

^2L{RmF) 


jiijk) 

3R(RmG) 


Aijk) 

•iL{RmIu) 


(dvr) 


m 


Rrr 


rrij 


Ml 




+ I 


m 


Rn 


mi 




Ml 


M] 


dH r d^q 


NX 


^ J {2xY J {2xY — Ml {i — qY — MI^^ — Mp 1“^ — m^^ 


M 


NXi 


dH f d^q i-q 


(11.17) 


(11.18) 

1 


AMpMp^ J (27r)4 J (27r)4 £2 _ mI {£ - qY - q2 _ mI £2 _ ^2 _ 


m 


2 ’ 


(11.19) 


M; 


d^l r d^q 


NX 


{2xYJ (27r)4 .^2 _ Ml g2 _ Ml^, £‘^ - £‘^ - {£ - qY - 


2 ’ 


( 11 . 20 ) 


M 


NXi 


dH r d‘^q £2 


MfMf, j (27r)4 J {2fY ^2 - Ml q^ - £‘^ - Ml £^ - {£ - qY - ml ’ 

( 11 . 21 ) 


with I{x) = X log x/(l — x). Note that in the derivation of the above formula, the CP phases except 
the Yukawa couplings are neglected. The contribution of the left diagram can be understood as 
linear seesaw like formula by splitting the diagram into two Dirac masses induced at one-loop level. 
For the center and right diagrams, there are two kinds of contributions coming from right and left 
chiralities of the internal fermions. In other words, these two contributions to the neutrino masses 
come from the masses or momenta of the F propagators in the loop respectively. The neutrino 
mass generation can be understood as follows. Due to the global U{1) symmetry breaking by the 
VEV of S, the mixing between rj, and x^' occurs. Then since the global 17(1) symmetry is 
correlated with the lepton number conservation, the 17(1) symmetry breaking implies breaking of 
the lepton number. Thus the neutrino Majorana mass term is generated after the 17(1) symmetry 
breaking. 

The neutrino mass matrix computed above can be diagonalized by the Pontecorvo-Maki- 
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FIG. 2: Numerical calculation of the loop functions where the other masses are fixed as Mp. = Mp^ = 1.5 
TeV and = mj^ = 1.2 TeV which are typical sample points to discuss DM and leptogenesis as we will 
see below. 


Nakagawa-Sakata matrix C/pMNS [76]; Upyi-j^^{mi,)UpMNS = diag(mi/^,mjyj,The neutrino 
masses and their mixing angles are measured by experiments [77], and these values depend on 
normal or inverted mass hierarchy. In our model, the order of magnitude of the neutrino masses 
can roughly be estimated as rriu ~ y'^Y'^O'^^^Iioopv' where Y is the dominant Yukawa coupling in 
UNx^ VNx'^ v'nx' y^Nx'^ Omix represents the mixing matrix of O^, , and Iioop is the loop function. 

Thus one can find that the order of ~ 10“® is required to obtain the experimental value 

niu ~ 0.1 eV with the typical assumed mixing angle Omix 10 /loop 0.1 and v' 10^ GeV. 
Note that the neutrino mass matrix should be proportional to the VEV v' since v' is the origin of 
the lepton number violation. 

We numerically compute the loop functions with the public code SecDec [78] in order to evaluate 
the neutrino masses more precisely in this model. Here one should note that the loop functions 
^ 3 R{R I ) ^3L{RaIb) ^ divergence. This is obvious from the definition of the loop 

functions in Eq. (11.20) and (11.21). However the divergent terms eventually cancel out with 
each other as follows. The loop functions can be regularized with dimensional regularization and 
expanded around dimension / = 4 which can be done within SecDec. With a brief evaluation, one 
can see that the terms including divergences are independent on at least either of m orn which is the 
index of the scalar mass eigenvalues. Thus the loop function including a divergence j ) can 

be written as j ) = )~^Taking into account this fact and the 

orthogonality of the mixing matrices and which means 0^0^^ = Ylm ^im^jm ~ 











one can see that the first two terms and vanish after taking the summation 

over n. Moreover the remaining third divergent term also cancels after all the relevant terms are 
summed in Eq. (11.15) and (11.16) as 
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Thus the divergent terms do not contribute to the neutrino masses. The numerical value of the 
loop functions are almost fixed by the maximum mass in Mp^, Mi^Xj, Mp^, '^RRm and and 
the result obtained by using SecDec is shown in Fig. 2.^ The numerical calculation shows that the 
loop function j ^ coming from the right diagram in Fig. 1 gives a dominant contribution to 

the neutrino masses. 


C. LFV processes 


We should take into account lepton flavor violations (LFVs) such as /r —?■ ey, which typically 
provide strong constraints on radiative neutrino mass models. In our case, such processes arise 
through only the term, and analyses are very similar with the case of the Ma model [79], and 
the Yukawa couplings y^x) VNx'i v'nx' y'Nx' constrained by the LFV processes at least at 

one-loop level. Among the LFV processes, we focus on the one-loop induced y —)• ey that gives the 
most stringent constraint on y^ and the mediating particles F and y"*". The resulting formula for 
y —7- ey and its experimental bound [80] are given by 

3 

Y.{y;)^l{yn)^2F2{^i) 


Br(y —)• ey) = 


ScTp 


647rG^m"^. 


i=l 


< 5.7 X 10 


-13 


(11.23) 


with F 2 {ii) = 


1 - 6 ^,+ 243+ 342-642111^, 


(11.24) 


6 ( 1 - 4 )^ 

where Gp = 1T7 x 10“® GeV“2 ig the Fermi constant, 4 = and aem = e^/{4:p) ~ 1/137 

is the electromagnetic fine structure constant. The simplest way to avoid this constraint is to 
assume y,, to be diagonal, since its formula is proportional to the off-diagonal elements of yr, as 
one can see in Eq. (11.23). In this case we expect the neutrino mixings can be derived through 
the other Yukawa couplings ywx) yiVx'; J/wx’ l^TVx" Otherwise the parameters are constrained as 
yrj < 0.01 and m^± > 200 GeV. For instance, with the values jjri = 0.01, m^± = 200 GeV and 
-^ 2 ( 4 ) = 1/6) tliG maximum branching ratio is found to be Br(y —)■ ey) ~ 1.4 x 10“33, 


At most 1% error is included in the numerical calculation. 
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Since we take < 0(0.01) for the LFV constraint and 0(1) TeV of the new particle masses 
in the loop, the other Yukawa couplings Vnx^ VNx'^ v'nx' y'Nx' t)e roughly larger than 10“^ 

in order to obtain the scale of the observed neutrino mass niy ~ 0.1 eV assuming Omix ~ 10 ^ as 
discussed in the previous section. 


D. Goldstone Boson 


Here we mention some issues on the GB. Due to the direct consequence of our global ?7(1) 
symmetry, the GB remains as a physical state, which could be constrained by some experiments. In 
our case, the constraint comes from the invisible decay of the SM Higgs boson, and its decay width 
can be computed with the coupling given in Eq. (H.9) to be r{h —)• GG) = sin^ Q;/(327ru'^). 
This decay width should be smaller than 1.2 MeV at 95% confidential level [81], and thus we get 
the constraint 


sina\ 

~Kr) 




< 2.5. 


(11.25) 


Moreover, sin a itself is constrained by the latest LHG searches by ATLAS and GMS to be (con¬ 
servatively) sin a < 0.2 [82]. Therefore the above constraint is translated to the constraint on the 
VEV as v' > 800 GeV. However since we take v' ~ 10^ GeV for successful leptogenesis, this bound 
is easily satisfied in our case. 

Another bound comes from the Supernova 1987A observations and simulations, which tell us 
the following relation [83]: 


lA$sl <0.011 


rriH 


500 MeV 


2 


(11.26) 


This bound also does not affect to our model seriously, since both mn and are taken to be 
free values of physical parameters. 


III. DARK MATTER 

We have two DM candidates which are the lightest fermion Fi and the lightest mass eigenstate 
of the scalars (ry, xoj Xo)r- These DM candidates can be stabilized by the accidental Z 2 symmetry 
but not a remnant symmetry of the global U{1) symmetry. This accidental Z 2 symmetry could be 
understood as a kind of the accidental symmetry which has been discussed for gauged U{1)b-l in 
Ref. [84]. 
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FIG. 3: Diagrams of the DM annihilations where DM is identified as the lightest scalar mass eigenstate. 


For the fermionic DM candidate Fi, the relevant coupling for DM annihilations is only the 
Yukawa coupling and the required strength of the Yukawa coupling is 0(1) for the DM mass 
above the electroweak scale in order to accommodate the observed DM relic density. On the other 
hand, small coupling < 0(0.01) is needed to evade the LFV constraint as has been discussed in 
the previous section. Thus the fermionic DM candidate Fi conflicts with the LFV constraint.® 
Therefore we identify the lightest mass eigenstate of the scalars as a DM candidate. The mixing 
angles among (ri,xo,Xo)R are induced by the scalar couplings A, A', A", and the magnitude of the 
mixing angles should roughly be Omix ~ 10“^ in order to reproduce the measured neutrino masses 
without conflict with the // —)■ ey process as discussed in the previous section. This order of the 
magnitude of the mixing angles can be achieved with the parameter setting given in Sec. II A. 
Since the full scalar potential given by Eq. (II.2) is rather complicated, we take into account only 
A, A', A", A$, A$s, As and A$^ for simplicity. Since the required order of the magnitude of the 
couplings A, A', A", A<i>s, As is very small, they would not affect to the computation of the DM relic 
density and detection probability. However the coupling A$s is important to induce the mixing 
angle sin a, and A$,j is responsible for direct detection of DM since this coupling generates the 
dominant contribution to the elastic scattering with nuclei mediated by the SM-like Higgs boson 

® Although the LFV constraint may be satisfied by considering a diagonal Yukawa matrix or specific flavor structure, 
we do not discuss this case. 
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FIG. 4: Numerical results in the (A$^, ttidm) plane for non-degenerate case (left plot) and 

degenerate case (right plot), where the parameters are fixed to be mn = 2 TeV and sin a = 0.1. 

h. 

The diagrams of the DM annihilations are shown in Fig. 3. The DM couplings in the scalar 
potential are basically weak in our parameter setting. However since the scalar DM candidate 
includes the SU{2)l doublet inert scalar r]ji, DM can annihilate into the gauge bosons via the 
gauge interactions in order to satisfy the observed DM relic density if the inert doublet scalar 
component of the DM candidate is sufficiently large. This can be achieved with a smaller (M^)ii 
compared to (M ^)22 and (M ^)33 in Eq. (II.6), and we consider such a case. In this case, the 
annihilation channels in the first line in Fig. 3 become dominant processes to determine the DM 
relic density. 

The DM relic density can be evaluated by using micrOMEGAs [85] and the results are shown 
in the (A$,y, m-DM) plane in Fig. 4 where the heavier CP-even Higgs boson mass is taken to be 
niH = 2 TeV and the mixing angle is sin a = 0.1 as an example. Here we take the negative 
values of the scalar couplings A, A', A". If the scalar couplings are positive, the lightest scalar in 
the imaginary components becomes DM candidate instead of the real components. The left plot 
shows the case that the masses of the other two heavier mass eigenstates are twice of the DM mass 
(lightest state), and the right plot shows the case that the heavier states are degenerate with the 
DM state. The red colored band represents the region satisfying the observed DM relic density 
by PLANCK within 3a confidence level [86]. The blue region is excluded by the direct detection 
experiment LUX [87], and the green region is expected to be tested by the future direct detection 
experiment XENONIT [88]. From these plots, one can see that when the lightest DM state is non¬ 
degenerate with the other heavier states, DM is close to the inert doublet DM (left plot), because 
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there is the mass threshold moM ~ 530 GeV for the left plot in Fig. 4, which is the same property 
of the inert doublet DM [89]. As well-known, the inert doublet scalar DM candidate can satisfy the 
observed relic density in the mass ranges of tudm ~ 60 GeV and ttidm ^ 530 GeV. On the other 
hand, the mass threshold can be lower as moM ~ 250 GeV if and X^' are degenerate with DM as 
one can see from the right plot in Fig. 4. This is because the interactions of the singlets and x^' 
are described by the scalar potential, and extremely limited in the case of the simplified potential. 
There is a small resonance feature at tudm ~ 1 TeV due to the channel DMDM —)• H* —>■ SMSM, 
however the resonance is not strong because of the small mixing angle sin a = 0.1. 


IV. RESONANT LEPTOGENESIS 


We consider the thermal leptogenesis in this model [90]. The lepton number asymmetry is 
expected to be generated through the out-of-equilibrium decay of the lightest Majorana fermions 
Ni and Xi, if we impose the lepton number of T as —1. Although the Yukawa coupling is 
required to be smaller than 0(0.01) from the LFV constraint, this is large enough that F and 
the SM leptons are in the thermal equilibrium. Thus, the generated lepton number asymmetry in 
the F sector can instantaneously be converted into the SM leptons, and then the baryon number 
asymmetry can be generated through sphaleron process. 

After the global U{1) symmetry breaking, the Yukawa interactions for Majorana fermions are 
written as 

C D F^{Ym,Pl + Y;,^Pr)N'x* + F^YxTl + Y^Tr)X'x* 

+ F^Yx^Pl + Y^Tr)N'x' + F^iYx^L + Y^^Pr)X’x' 

+ F(V4Pl + Yx,Pr)N'x^* + F{Y^Tl + Yx,PR)Xix'* 

+ F{YkPL + Yn.Pr)N'x + F{Y'x^Pl + Yx,Pr)X'x, 

where N- and Y( are expressed as the mass eigenstates of each Majorana fermion. Hereafter, we 
abbreviate them to Ni and Xi for convenience. The Yukawa couplings are redefined as 


YNix = VNix^osOi, Y^.^ = y'jx.^smOi, 
YXiX = = y'xix^osOi, 

YXix' = VNiX' ^kx' = ykx' 

Yxix' = y*Xx' ^kx' = ~ykx' 

where 9i is the mixing angle of the i-th generation. 


(IV.2) 
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FIG. 5: Lepton number violating decay and scattering processes for leptogenesis where A and A' represent 
the new scalar singlets, or x^' and li {i = 1 — 2) is Ni or Xi 


In this model, we consider that the TeV scale masses of Majorana fermions and v' r\j 10^ GeV 
so that the annihilation channels A^iiVi, AiXi —)• GG are decoupled from the thermal bath before 
the temperature of the universe T ~ 10 TeV. Otherwise the generated lepton asymmetry would 
be washed out by these lepton number violating processes. 

Although the TeV scale of the masses seems to be too small to realize the sufficient CP asym¬ 
metry associated with the decay processes for the generation of the required baryon number asym¬ 
metry, the generated baryon asymmetry can be enhanced by the resonance effect as known in 
resonant leptogenesis [91-99]. We define the parameter e as the amplitude of the CP asymmetry. 
The dominant contribution for the CP asymmetry comes from the interference between the tree 
diagram and the one-loop self-energy diagram as depicted in the upper line in Fig. 5 and its formula 
is given by 


{Ml - Ml)Mj,Ti, 
{Ml-Mlf + MlT]J 


(IV.3) 


where Mj. and Tj. are the mass and the decay rate of (/ = V or V) respectively. From this 
equation, the maximum enhancement is caused when Ml — Ml = Mj^Tj^. In this model, we can 
take a larger T compared to that of the second lightest right-handed neutrino in the canonical 
resonant leptogenesis at TeV scale, since the Yukawa couplings can be large without conflicting 
with the observed neutrino masses due to the loop suppression. Thus, the required magnitude of 
the degeneracy of the Majorana fermion mass can be quite milder to generate the sufficient baryon 
number asymmetry than those in the usual resonant leptogenesis at TeV scale [73, 100]. One may 
think that the baryon asymmetry should be correlated with the VEV of the singlet scalar S since 
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the B — L breaking occurs with only v'. Indeed this B — L breaking effect is included in the 
total mass matrix of Ni, Xi and the active neutrinos t'j. For example, the B — L violating Dirac 
mass term between r'j and Xj is induced at one-loop level as can be seen from the left diagram 
in Fig. 1. Thus the effect of the B — L breaking is included in the masses of Ni and Xi, and one 
can understand that the baryon asymmetry is generated through the breaking effect in the mass 
matrix. 

We need to take into account washout effects to evaluate the baryon number asymmetry. The 
generated lepton number asymmetry could be washed out through the lepton number violating 
2-2 scattering processes and the inverse decay of /*. However, if the relevant Yukawa couplings are 
small enough, these processes can be nearly decoupled before the temperature of the thermal plasma 
decreases to T < Mj^. Thus, the washout of the generated lepton number asymmetry is expected 
to be suppressed sufficiently in this period. In order to examine this quantitatively, we numerically 
solve the coupled Boltzmann equations for the number density of A^i, Xi and the lepton number 
asymmetry. We introduce the number density of Ni, Xi and the lepton number asymmetry in 
the comoving volume as Y^i = Yxi = and Yp = {np — np)/s respectively, by using 

the entropy density s and the number densities which are expressed by nxi, nx^, np and np. 
Their equilibrium values are given by Y^^ = z^K 2 {z), where z is defined hy z = 

is the number of relativistic degrees of freedom and K 2 {z) is the modified Bessel function of the 
second kind with the order 2. Since we assume Yp is immediately translated into the SM leptons 
as we mentioned above, we use the relation B = ^{B — L) which is derived from the chemical 
equilibrium condition in this model, and the baryon number asymmetry Yp in the present Universe 
is estimated as Yp = —§^Yp[zp-sx)^ where zpxi = -^/i/Few is related to the sphaleron decoupling 
temperature Tew- 

The coupled Boltzmann equations for the leptogenesis in our model are written as [101] 


dYx, 


Yn, 




dz sH \Y^l 


_ 

dz sH I 


dYx, 

dz 

dYp 

dz 


-1 hx, 


z 

2Yp 

~YP 


ex 


Y^ 


- 1 I 7^1 + ex 




-1 hw 


E 


iNj + IXj 

4 


+ ^ ilFAFA' + Iff A A') 
A,A' 


(IV.4) 

(IV.5) 


(IV.6) 
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where is defined by 


* TT^ 


K^ 


Mu 

~Y 


pD 


(IV.7) 


with the modified Bessel function of the second kind Ki[z) with the order 1, UaUj is the reaction 
density for the scattering process ab uu ij which is given by 


T f°° (y/s\ 

labij — / ds d'abijis)y/sKi f j , (IV.8) 

with Smin = max [{nia + rnb)"^, (m* + mj)^] and the reduced cross section dabij{s). There are two 
kinds of the processes TA -H- FIS' and FF uu A A' for the scattering processes as depicted in 
the bottom line of Fig. 5. The reduced cross section in our model is rather complicated since 
a lot of particles exist, but can be straightforwardly computed from the Lagrangian as same as 
Ref. [102, 103]. We solve the coupled Boltzmann equations numerically. 

The decay of the lightest Majorana fermions should be out of thermal equilibrium so that 
the lepton number asymmetry can be generated through their decays. If we express the Hubble 
parameter as H, this condition is given by H > Tp at T ~ Mj^. Since we assume that each of 
the Majorana fermion mass is = 5 TeV, Mj^ = Mp(l + 6M) and M/g = 6 TeV, the Yukawa 
couplings of Ni and Xi should be 0(10“®). On the other hand, the rest of the Yukawa couplings 
should be 0(10“^) in order to generate the appropriate neutrino masses. Here we set Yukawa 
couplings in Eq. (IV.2) to be yux = 1-5 x 10“® and y/j.sx ~ 10“^ and 9i = We show the result 
for 6M = 10“® in Fig. 6 as an example, and also the generated baryon number asymmetry for each 
value of 5M and Mq in Fig. 7. Through this analysis, the masses of F, x and x' are fixed to be 
Mp = 1.5 TeV and My, = M^/ = 1.2 TeV respectively. These parameter set satisfies the condition 
for the DM phenomenology we discussed in the previous section. 

From the left panel in Fig. 6, we can see that the required baryon number asymmetry Yp can 
be obtained in this model. In the right panel, we plot the behavior of the relevant reaction rate 
for each process. This panel shows that the reaction rates of the inverse decay process and the 
lepton number-violating process induced by the s-channel exchange are quite large for a long 
time. Thus the baryon number asymmetry cannot be generated quickly until rather a late period. 
After T ~ M/^, the generated baryon number asymmetry gradually increases and then the required 
value can be realized. This is because these processes are suppressed by the Boltzmann factor. 

We show the relation between the generated baryon number asymmetry and the mass degen¬ 
eracy of Majorana fermions in Fig. 7. Notice here that the generated baryon number asymmetry 
is always smaller than the required value in the case ~ 3 TeV. In this model, we can realize 
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FIG. 6: Left panel: Evolution of Yb, and Yxi where the horizontal black dashed line represents the 
required value of the baryon number asymmetry. Right panel: Reaction rates T/H oi the processes that 
have crucial effects for the baryon number asymmetry. 

10 '^ 

10 -^° 

10 - 2 ° 
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M, =4 TeV X 
M||=5 TeV X 


FIG. 7: 6M dependence of the generated baryon number asymmetry Yb with the different value of M/^. 
The value of Mj^ is fixed to 3 TeV, 4 TeV and 5 TeV respectively. 

the large Yukawa couplings to explain the small neutrino mass due to the two-loop effects and 
then T becomes larger compared to tree and one-loop neutrino mass models. Thus, the required 
mass degeneracy can be milder. However, the large Yukawa couplings cause the large washout 
effects. Since the Boltzmann suppression does not work well in the case of small the large 
washout effects remain until quite a late period compared to the heavier cases. Thus the most of 
the generated baryon number asymmetry is washed out. We find that the observed baryon number 
asymmetry can be generated when the mass degeneracy is 5M = (10“® — 10“^) for Mj^ =4 — 5 
TeV, as can be seen in Fig. 7. As we mentioned above, the magnitude of this mass degeneracy is 
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quite milder than the canonical seesaw case for each value of due to the loop effects.® 

V. CONCLUSIONS 

We have studied a two-loop induced radiative neutrino model at TeV scale with a C/(l) global 
symmetry, in which two types of DM candidates (the lightest one of fermion or scalar) can be 
involved. The loop-induced neutrino masses have been evaluated appropriately and phenomenology 
of DM has also been discussed. The fermionic DM candidate is disfavored if we reproduce the 
measured neutrino masses and take into account the constraint from LFV with the same order of 
all the elements of Then we have found that the scalar can be a good DM candidate which 
satisfies the observed relic density and the DM direct detection bound. We also found that the 
direct detection rate of DM is controlled by the coupling and some parameter region can be 
testable by the next future direct detection experiment XENONIT. 

We have discussed baryon number asymmetry through the resonant leptogenesis with multi¬ 
sources scenario, in which the large Yukawa couplings (that is required to compensate the tiny 
neutrino masses at two-loop level) make the large CP asymmetry, but also cause the large washout 
effects. We have shown that the required baryon number asymmetry can be obtained for the 
parameter region i.e., 5M = (10“^ — 10“^) for =4 — 5 TeV, where I = N or X. In this case, 
the lightest Majorana fermions should satisfy M/j > 3 TeV to suppress the washout effects by the 
Boltzmann factor. For larger , the required magnitude of the mass degeneracy is rather milder 
than the canonical seesaw case even at TeV scale models. 
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